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1. Introduction 

Recent tests of the unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) matrix have indi- 
cated some tensions at the 2-3a level [|l]-|h. In many cases, the constraints on CKM unitarity are 
limited by the precision of the theoretical inputs, in particular the hadronic matrix elements that 
characterise the strong interaction dynamics of weak processes. It is therefore imperative that these 
matrix elements are determined as precisely as possible. 

The HPQCD collaboration is currently undertaking a suite of precision lattice calculations of 
heavy-light mesons to reduce the uncertainties associated with the theoretical inputs into CKM uni- 
tarity fits. New calculations of the decay constants and /g ( using the Highly Improved Staggered 
Quark (HISQ) action reached a precision at the 2% level [|], §]. These results are the most precise 
currently available for these decay constants and were made possible by the chiral properties and 
reduced taste-breaking uncertainties of the HISQ action. 

In [^J, the heavy-light currents were computed using HISQ light quarks and nonrelativistic 
QCD (NRQCD) heavy quarks. This calculation requires matching the heavy-light axial-vector and 
vector currents in the effective theory on the lattice with full QCD. In this article we report on the 
one loop perturbative matching of the HISQ-NRQCD currents for massless HISQ quarks. We also 
discuss the extension to heavy-heavy currents. 

In the next section we describe the quark and gluon actions used in our calculation. We then 
review the matching formalism for heavy-light currents. In Section |] we present our results for a 
range of heavy quark masses and include the quark renormalisation parameters. We discuss the 
extension to heavy-heavy currents in Section || and conclude in Section ^. 



2. Lattice Actions 



We use the Symanzik improved gluon action with tree level coefficients [0], given by 



S G 



JU Q x,jl>V 



4u Q 



(2.1) 



Here P^ v is the usual plaquette and R^ v the six-link loop, with j8 = 2N c /g 2 and uq a tadpole 
improvement factor [0]. Radiative improvements to the gluon action do not contribute to our one 
loop matching calculation, because one loop radiative improvement generates an 0{oc s ) insertion 
in the gluon propagator that would only contribute at higher orders. 

We include a gauge-fixing term and, where possible, we confirm that gauge invariant quantities 
are independent of our choice of gauge parameter by working in both Feynman and Landau gauges, 
which we denote t, = 1 and E, = respectively. 

We use the Highly Improved Staggered Quark (HISQ) action [Jsj] for the light quarks. The 
HISQ action has been used successfully to simulate both b and c quark systems [^, 10]. Taste- 
breaking discretisation errors are significantly reduced by two levels of link fattening and ff(a 4 m.g) 
uncertainties are suppressed by powers of v/c (where v is the quark velocity) through a tuned a 
coefficient for the three-link "Naik" term flsll . We write the action as 



5HisQ = « 4 I?W(7^V» ISQ + mo)vAW, with V» IS< 3 = V 



(FUF) 



(2.2) 
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The superscripts indicate that the first operator, Vif UF \ is built from the full HISQ-smeared links: 

^ = (^-Z^r)*^ where ^ = nf 1 +^) (23) 

V P+V / P+V V / symmetrised 

and ^ is a reunitarisation operator. The second operator in Equation (|2.2|), vjf 7 ', includes only 



one level of smearing. We work with massless quarks, so we set the bare quark mass, mo, and the 
tuning parameter, e, to zero. 

For the heavy quark fields, x ¥(x,t), we use the NRQCD action of JTTj] , which is improved 
through /Mq) and G(a r) and includes the leading relativistic <^(1/Mq) correction. The full 
NRQCD action is 

t \( aSH\ ( aH \" t / aH Q \ n / aSH\ 

where 1/// = y/^(x,f) and = y/(x,f — 1). The leading kinetic term is given by 

A <2) 

"»=-2^o' (25) 
and 5// includes higher order improvement terms, full details of which are given in, for example, 



[11]. We use the tree level values of c, = 1 for all the coefficients, Cj, of the higher order operators 



in 8H, and do not consider the effects of radiative improvement of the NRQCD action. 

3. Matching Procedure 

On the lattice, the heavy-light axial-vector and vector current operators mix with higher order 
operators under renormalisation. We relate the lattice and continuum currents perturbatively and 
extract the mixing matrix elements at one loop. For massless HISQ quarks, the results for axial- 
vector and vector currents are identical. Our strategy for the perturbative matching of heavy-light 
currents follows that developed in [ 12] and outlined in [p"3|]. 



We require three lattice currents to match the temporal component of the vector current through 

&(a s ,AQCD/M Q ,a s /(aMo),a s AQcu/M ). These are 

= q(x)T^Q(x) , = -^-L-^(x)r M7 . ^Q{x) , J® = --J—^( x ) Y .^ Y(i r^Q(x). 

(3.1) 

Here the Q fields are four component Dirac spinors with the upper two components given by the 
two component NRQCD field and lower components equal to zero. The operator represents the 
vector current operator, T^ = y il . 

The matrix elements of the vector current in full QCD are related to those in the effective 
theory via 

(V ) = (l + c, Po W ) (Jf > + (l + a,*™) (J^ b ) + (l + c,pf 0) ) (4 2) ' SUb ); (3-2) 
(^)=fl + «.vPo ra )(4 0) ) + (4 1) ' SUb )- (3-3) 
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Here we have expressed the lattice currents in terms of the subtracted currents, 

J^^jV-a^, (3.4) 
for i = 1, 2. The subtracted currents are more physical and have improved power law behaviour 



[12]. Note that, for the spatial components, we match through G((x s ,Aqct>/Mq,(X s / (oMq)). 
We match at zero external quark momentum; the matching coefficients are given by 

- l ^j-\{C q + C H )-C\ (3.5) 

-l£) -\{C q + C H )-C M -^ ) -C\ (3-6) 
4 Vo) , (3-7) 
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(3.8) 



where the contributions from continuum QCD are given in [12, 13]. The renormalisation parame- 
ters C q , Ce and Cm are the one loop HISQ wavefunction renormalisation and the NRQCD wave- 
function and mass renormalisation parameters respectively. We have written the pole mass, which 
is common to both lattice and continuum theories, in terms of the bare quark mass and must there- 
fore include the one loop mass renormalisation in p\. 



The Qj in Equations (3^) to (^8) are the one loop mixing matrix elements that arise from 



the mixing of the currents. The matrix element £ 2 includes a term that removes an 0(aa s ) discreti- 



sation error from Jfr' [ |T2[ 13[ ]. Thus our matching procedure ensures that ff{a s /Mo) and G[aa 



corrections are made at the same time. There is a second dimension four current operator that is 

(2) H — I I — h 

equivalent to Jq via the equations of motion [12, 13], which we include in the determination of 
Ca (for/ = 1,2). 

4. Results 

We calculate the mixing matrix elements and renormalisation parameters with two indepen- 
dent methods. Our first method used the automated lattice perturbation theory routines HiPPy and 



HP src [14]. We performed these automated lattice perturbation theory calculations on the Darwin 
cluster at the Cambridge High Performance Computing Service and the Sporades cluster at the 
College of William and Mary with routines adapted for parallel computers using MPI (Message 
Passing Interface). Our second method used Mathematica and Fortran routines to construct 
the appropriate Feynman integrands; these were then evaluated with Vegas JIB]]. 



We tested our code in a number of ways. We reproduced the NRQCD-AsqTad results of [ |13[ ] 
and in many cases, we established that gauge invariant quantities, such as the mass renormalisation, 
were gauge parameter independent by working in both Feynman and Landau gauges. We also 
confirmed that our results exhibited the correct infrared behaviour, regulating any divergences with 
a gluon mass and using subtraction functions to ensure divergences were correctly handled by 
Vegas . We believe that our two methods were sufficiently independent that agreement between 
these methods provides a stringent check of our results. 
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nM r r r>( y °) o^o) n (y ) r (v ) (V k ) r {v k ) 



2.688 -0.360(2) 0.262 -0.108(2) 0.012(2) -0.712(4) -0.1144 -0.034(2) 0.0382 

2.650 -0.370(2) 0.267 -0.111(2) 0.013(2) -0.693(4) -0.1157 -0.034(2) 0.0386 

1.832 -0.658(2) 0.466 -0.162(2) -0.042(3) -0.314(4) -0.1593 0.018(2) 0.0532 

1.826 -0.662(2) 0.468 -0.164(2) -0.043(3) -0.311(4) -0.1595 0.020(2) 0.0532 

Table 1: One-loop renormalisation parameters for the heavy-light vector current in Feynman gauge. For the 
p/ V,u ' the quoted uncertainties are the errors from each contribution added in quadrature, whilst for Ch, Cm 

and the Cio^ the uncertainty is the statistical error from numerical integration. Unless otherwise indicated, 
the uncertainties are 1 or smaller in the final digit. 

Many of the one loop parameters that we calculate are infrared divergent. We decompose our 
results into into an infrared finite term and an infrared divergent contribution, which we denote 
with a superscript IR . Thus we write 

Z q = 1 + a, (Cf + C q ) + G(a]) , Cf = -L [l + - l)] log (a 2 A 2 ) , 

Z H = 1 + a s (Cf + C H ) + G(a]) , Cf = -L [-2 + (§ - 1)] log (a 2 X 2 ) , (4.1) 

Coo = Coo + + ^(«.v ) , C = ^ log(a 2 A 2 ). (4.2) 

Here A is the gluon mass, introduced to regulate the infrared behaviour, and <^ is the gauge param- 
eter. We confirm that both the gluon mass dependence and all infrared divergences cancel between 
the lattice and continuum one loop coefficients. Note that the NRQCD mass renormalisation, 
Zm = 1 + a s C M + ^(« 2 ), is infrared finite. 

We tabulate results for the infrared finite contributions to the renormalisation parameters, mix- 
ing matrix elements and matching parameters for the heavy-light vector current at four different 
heavy quark masses in Table ||. For the NRQCD action we present results with c,- = 1 and sta- 
bility parameter n = 4. We use the Landau link definition of the tadpole improvement factor, 

= 0.7503(1). Only the matching coefficient p\ has a tadpole correction coefficient. This 
correction contributes to £^ ' and is given by = Uq\ The HISQ wavefunction renormali- 
sation is of course independent of the NRQCD mass and we find C q (t, = 0) = 0.1145(1) and 
Q(§ = 1) = -0.3940(1). 

5. Heavy-heavy Currents 

We are currently extending our calculation to heavy-heavy currents, with non-zero HISQ mass. 
Moving from massless to massive relativistic quarks complicates the matching procedure. In the 
former case, quarks and antiquaries at zero spatial momentum are indistinguishable and conse- 
quently scattering and annihilation processes give identical results. In the massive case, however, 
we must distinguish between quarks and antiquaries. Massive HISQ quarks also complicate the nu- 
merical integration considerably. The chief difficulty arises for the Aq and annihilation currents, 
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which contain a Coulomb singularity that must be handled with a subtraction function. Moreover, 
in the automated perturbation theory routines, the pole in NRQCD propagator crosses the integra- 
tion contour and we must introduce a triple contour to ensure the stability of numerical integration 

ph. 



6. Summary 

We have determined the one loop matching coefficients and renormalisation parameters re- 
quired to match the axial-vector and vector currents on the lattice to full QCD. We used the mass- 
less HISQ action for the light quarks and NRQCD for the heavy quarks and match at zero external 
quark momentum. 

These matching coefficients are important ingredients in the determination of heavy-light 
mesonic decays in lattice QCD studies Studies of the B s meson using the relativistic HISQ 
action for both b and s quarks have recently been carried out [|J]. Such an approach has the advan- 
tage that perturbative matching, which is generally the dominant source of error in the extraction 
of decay constants, is not required. Computations at the physical b quark mass are currently pro- 
hibitively expensive, however, and require an extrapolation up to the b quark mass. Furthermore, 
simulations of the B meson are not yet feasible, because the use of light valence quarks and close- 
to-physical b quark masses require large lattices with fine lattice spacings. Thus an effective theory 
approach to heavy-light systems remains the most efficient framework for precise predictions of 
IbJ/b and fs- Such calculations require the perturbative matching calculation reported in this 
article. 

The extension to heavy-heavy currents will enable the extraction of form factors for the B — > 
D^iv semileptonic decays and the leptonic decay of the B c meson from lattice NRQCD compu- 
tations. These calculations are underway. 
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